Cross-link governed dynamics of biopolymer networks 
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Networks of stiff biopolymers cross-linked by transient linker proteins exhibit complex stress 
relaxation, enabling network flow at long times. We present a model for the dynamics controlled by 
cross-links in such networks. We show that a single microscopic timescale for cross-linker unbinding 
leads to a broad spectrum of macroscopic relaxation times and a shear modulus G ~ u 1 ^ 2 for 
low frequencies uj. This model quantitatively describes the measured rheology of actin networks 
cross-linked with a-Actinin-4 over more than four decades in frequency. 



PACS numbers: 87.16.Ka,83.80.Lz,87.15.La,87.15.H- 

Reconstituted biopolymers such as actin are excellent 
models for semi- flexible polymers, with network mechan- 
ics and dynamics that are strikingly different from flex- 
ible polymer networks [THE]. One essential feature set- 
ting biopolymer networks apart from rubber-like mate- 
rials is the intrinsic dynamics of their cross-links. Such 
systems represent a distinct class of polymeric materials 
whose long-time dynamics are not governed by viscosity 
or reptation [9], but rather, by the transient nature of 
their cross-links. This can give rise to a complex me- 
chanical response, particularly at long times, where the 
network is expected to flow. Such flow can have impor- 
tant implications for cells, where their internal networks 
are constantly remodeling, reflecting the transient nature 
of their cross-links [TO] . The simplest possible descrip- 
tion of a material that is elastic on short timescales while 
flowing on long timescales is that of a Maxwell fluid; 
this exhibits a single relaxation time r, as depicted in 
Fig.[T] Indeed, some recent experiments on transient net- 
works have suggested the existence of a single relaxation 
time [11 ; by contrast, other experiments — probing longer 
time-scales — evince a more complex viscoelastic behav- 
ior [12J[T3], indicative of multiple relaxation times. Thus, 
the basic physical principles governing such transient net- 
works remain a mystery. A predictive theoretical model is 
essential to elucidate the effect of dynamic cross-linking, 
and to help explain the complex viscoelastic behavior ob- 
served experimentally. 

Here, we develop a microscopic model for long-time net- 
work relaxation that is controlled by cross-link dynam- 
ics. This cross-link governed dynamics (CGD) model de- 
scribes the structural relaxation that results from many 
independent unbinding and rebinding events. Using a 
combination of Monte Carlo simulations and an analytic 
approach, we demonstrate that this type of cross-link dy- 
namics yields power-law network rheology arising from a 
broad spectrum of relaxation rates. Our predictions are 
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FIG. 1: (Color online) A schematic of the frequency depen- 
dent shear modulus G* = G' + iG" . Non-permant networks 
can exhibit a response ranging from a single timescale (r) 
Maxwell-like behavior (blue lines) to a powerlaw regime with 
an exponent < 1 governed by a broad distribution of relax- 
ation times (> r) (red lines). Upper inset: for times longer 
than the unbinding time r Q ff , large scale conformational relax- 
ation can occur via linker unbinding (open cirlce) and subse- 
quent rebinding at a new location. Lower inset: for shorter 
times, only small-scale bend fluctuations between cross-links 
can relax, resulting in a plateau in G' for frequencies > l/r ff. 

in excellent quantitative agreement with experiments per- 
formed on actin networks cross-linked with the transient 
linker protein a-Actinin-4. 

The CGD model can be qualitatively understood in 
simple physical terms. We assume each filament to be 
cross-linked to the network, with an average spacing £ c . 
Only filament bending modes between cross-links can re- 
lax (Fig. [I] lower inset), and the thermalization of these 
modes results in an entropic, spring-like response. To 
account for transient cross-linking, we assume that the 
linkers unbind at a rate l/r Q ff (Fig. [I] upper inset); this 
initiates the relaxation of long wavelength (> £ c ) modes, 
giving rise to a softer spring constant and reduced macro- 
scopic modulus. However, the relaxation of successively 
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longer wavelength modes becomes slower, as an increas- 
ing number of unbinding events are needed for such a 
relaxation. This simple physical picture suggests a broad 
spectrum of relaxation times, as opposed to the single re- 
laxation time of the Maxwell model. As outlined below 
both simulations and an analytic treatment of this model 
yield power-law behavior with G ~ uu 1 / 2 below the char- 
acteristic frequency ujo — 2tt/t ^ (Figs. [2] A and B). 

We compare the basic predictions of this model to the 
rheology of a representative transiently cross-linked actin 
network. As a cross-linker, we use a-Actinin-4 [T3l fl4] . 
whose unbinding time r Q ff is reported to be in the range 
1 — 10s, similar to that of other biological cross-linkers. 
These gels [I2j [15] exhibit a low-frequency elastic shear 
modulus G' with a pronounced decay over three decades 
in frequency, while the viscous modulus G" exhibits a 
broad local maximum located near the characteristic fre- 
quency of cross-link unbinding [TTJ [12j [15] (Fig. [2^>). In 
the asymptotic low- frequency range, both moduli exhibit 
power-law rheology with an approximate exponent of 1/2, 
in agreement with our predictions. Such behavior clearly 
indicates a more complex stress relaxation than captured 
by the Maxwell model, which is governed by a single re- 
laxation time (Fig. [I]). Taken together, the theoretical 
and experimental results demonstrate a distinct cross-link 
governed regime of network dynamics. 

To develop a predictive microscopic model, we first con- 
sider a single polymer within the network, and then ex- 
tend the description to the macroscopic level. On length- 
scales longer than i c , the motion of the polymer is con- 
strained by its cross-linking to the surrounding network 
(see insets Fig. [I]). When a linker unbinds, a local con- 
straint is released, allowing for the thermal relaxation of 
the freed segment. This thermal relaxation occurs within 
a time r eq , which is typically of order milliseconds [3][4j[7]. 
We assume that this process is completed before the seg- 
ment rebinds to the network at a new location; thus, 
r eq <C r on , where r on is the rebinding time of the linkers. 
Furthermore, assuming r on <C r Q ff, only a small fraction 
of cross-links will be unbound at any given time, and si- 
multaneous unbinding of neighboring cross-links can be 
neglected. This suggests a coarse-grained description on 
length-scales longer than i c , in which independent un- 
binding events occur at a rate l/r ff . Since the relaxation 
of wavelengths less than £ c occurs at a much faster rate 
l/r eq , we use the worm-like chain model, where the equili- 
brated short wavelength fluctuations manifest themselves 
as an entropic stretch modulus /i t h ~ n 2 /£^.k^T [2j [4]. 
Here, n is the bending rigidity, &b is the Boltzmann con- 
stant and T is the temperature. In this description the 
coarse-grained energy is given by 

^ = }£[f |At n | 2 + ^(|Ar„|-4) 2 ], (1) 

c n 

where the sum extends over all cross-link positions r n , t n 
is the unit tangent vector, and A represents a discrete 
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FIG. 2: (Color online) A) The simulated rheology for fre- 
quencies below cjo- The shear modulus is normalized by the 
elastic plateau value Go . The inset shows the total power- 
spectrum C(uj) (blue circles) of distance fluctuations, as well 
as the fraction Cm coming from effective stretch fluctuations 
originating in undulations on length scales shorter than the 
cross-linking distance. The distance fluctuations are deter- 
mined over a length 16^ c of a polymer with a persistence length 
£ p — 32£ c and a total length 32£ c . The solid red line represents 
our analytical mean- field prediction. B) Measured linear rhe- 
ology of a 23.8 /xM actin network cross-linked with 0.238 /iM 
a-Actinin-4. The low frequency behavior is consistent with 
G rsj (zcj) 1//2 . The solid and dashed lines are global fits utiliz- 
ing our mean-field CGD model for the low frequency regime 
together with the known high frequency response [3j [4] . 

difference. 

Using Hcg, we study the dynamics arising from multi- 
ple linker unbinding and rebinding events, by performing 
2D simulations of a single polymer. An initial chain con- 
formation with periodic boundary conditions is randomly 
drawn from a Boltzmann distribution. Cross-link unbind- 
ing events are independent and result in the complete 
thermal equilibration of the two neighboring polymer seg- 
ments. This is numerically implemented via a Metropolis 
Monte Carlo algorithm. 

These simulations allow us to determine the equilibrium 
fluctuations of a single polymer, treating its surround- 
ing network as a rigid medium. According to the fluc- 
tuation dissipation theorem (FDT), the linear mechan- 
ical response of the polymer is encoded in the fluctua- 
tions of the extension, 5£, of the polymer. Interestingly, 
the simulations demonstrate that the power spectrum 
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C(u) = (\S£(uj)\ 2 ) depends on frequency as a fractional 
power-law, as shown in the inset of Fig. [2]A., indicating 
a broad underlying distribution of relaxation times. The 
exponent is consistent with —3/2 over five decades in fre- 
quency. Although this exponent also arises in the Rouse 
model for flexible polymers due to the viscous dynamics 
of longitudinal stretch modes [9], this is not the origin of 
the behavior found here. Our model does exhibit effec- 
tive longitudinal stretch modes; however, their contribu- 
tion Cm to the fluctuation spectrum is subdominant, as 
shown in the inset of Fig. [2}^. This demonstrates that the 
polymer's response to an applied tension is dominated by 
the dynamics of long wavelength transverse modes. 

The dynamical description of a single polymer can be 
extended to the network level by assuming that the net- 
work deforms afflnely. The macroscopic shear modulus 
G* is then related to the complex response function x 
of relative length extension of a single polymer in re- 
sponse to a tensile force: G* = p/(15x), where p is 
the length of polymer per unit volume [3j [4]. Ignoring 
end effects, the relative extension 5£/£ of a polymer seg- 
ment of length £ is conjugate to the uniform tension /, 
with S£(uj)/£ = xM/M- We use the FDT to calculate 
the imaginary part of the extensional response function 
£x"(oo) = uu(5£ 2 (uo)) /2k^T . Using a Kramers-Kronig re- 
lation, we compute the full complex response function x, 
required to obtain the network shear modulus [21 . Re- 
markably, below cj ff, the shear modulus depends on fre- 
quency as a power-law with an exponent of 1/2 (Fig. ^ 
consistent with experiments (Fig. [2^). 

To obtain further insight into this behavior, we develop 
a continuum analytical treatment of this model. We cal- 
culate the polymer displacement due to the unbinding and 
subsequent rebinding of a linker to the n-th cross-link site 
along the chain backbone. To capture this equilibration 
step, we approximate the thermal distribution as Gaus- 
sian and centered around the mechanical equilibrium of 
the coarse-grained chain. This allows us to separate the 
equilibration step into a deterministic move to the mini- 
mum energy position together with a stochastic thermal 
contribution. The deterministic part of the displacement 
5r n = — r$ from the initial (i) position to the local 
equilibrium position (eq) is determined by 
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This condition replaces the usual equation of motion bal- 
ancing drag and conservative forces in the low-Reynolds 
number regime. In this continuum long- wavelength de- 
scription, the leading order evolution equations are 
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where and ry are the transverse and longitudinal de- 
flections of the polymer with respect to its average direc- 
tion x. The noise £j_ captures thermal effects, including 
local bucking due to thermally-induced compression [16]. 
It can be calculated from a quadratic expansion of Hcg 
around its local mechanical equilibrium. For an inexten- 
sible polymer, the longitudinal component of this noise is 
subdominant and has been neglected. 

Importantly, the noise £j_ depends nonlinearly on the 
local state of the polymer and couples Eqs. (3|4). To 
explore this coupling, we artificially reduce the stretch 
modulus fi. In the limit \i <C /ith, the equations decou- 
ple and become exactly solvable, and the resulting trans- 
verse contribution to the fluctuation spectrum approaches 
C± ~ cj -7 / 4 . This can also be seen in our simulations 
with variable \i < /i t h in Figs. [3] A and B. As \i is re- 
duced below /i t h, Cj_ evolves toward C± ~ co> -7 / 4 , which 
can be seen by the flattening of the normalized spectrum 
in Fig. [3j3. In the limit /i <C /ith, the transverse bending 
dynamics are effectively those of a stiff filament fluctu- 
ating in a viscous solvent, for which the time-dependent 
fluctuations are (\S£(t)\ 2 ) ~ t 3 / 4 [S H [17]. Only in this 
decoupled limit, can one understand the dynamics within 
the framework of an effective viscosity provided by the 
transient cross-links [T6] . 

The nonlinear nature of the effective noise precludes 
a full analytical solution of the model in the limit of an 
inextensible polymer. Instead, further insight is gained by 
approximating the amplitude of the effective noise term 
by its mean-field value calculated from the equilibrium 
fluctuations of the polymer. In this approximation, the 
stochastic contributions are uncorrelated in both time and 
space, resulting in the response function [T6] 
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This response function captures the cross-link governed 
dynamics dominating on timescales > r ff. Further, we 
calculate the mean-field correlator, Cmf ~ co> -3 / 2 , in good 
agreement with the simulations presented in the inset of 
Fig. [2] As a further test, we perform simulations over 
a wide range of k and polymer lengths L; the predicted 
amplitudes are in excellent agreement with simulated am- 
plitudes, as shown in Fig. [3p. This further validates the 
assumptions made in our analytical approach. 

To obtain a description of the mechanical behavior on 
all time-scales, we combine the response due to cross-link 
governed dynamics with drag-limited short-time dynam- 
ics [3[ H]. Remarkably, the model can be fit to the ex- 
perimental data — over the full range of frequencies — with 
three parameters: the plateau modulus Go = 53Pa, the 
equilibration rate r eq = 0.02s and the unbinding time 
Toff = 3.2s (see Fig. [2]). This provides strong evidence 
that the low-frequency rheology of actin networks with 
the physiological linker a-Actinin-4 is governed by the 
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FIG. 3: (Color online) The power spectrum C±(uj) of lon- 
gitudinal fluctuations originating from transverse undulations 
for length scales longer than £ c , multiplied with u 3 ^ 2 (A) and 
cj 7 / 4 (B) for a range of polymer backbone compliances. C) The 
simulated amplitude of the power spectrum C(u))/r£c plotted 
against the 2D mean-field prediction for a range of polymer 
lengths and bending rigidities. 

linker-controlled dynamics. Furthermore, the fitting pro- 
cedure yields r eq <C r fr; this, together with the quality of 
the fit, lends credence to the separation of timescales as- 
sumed in our model. Such a separation of timescales also 
implies that the fluid viscosity does not affect the dynam- 
ics or rheology in the linker-governed regime, consistent 
with observations in other experiments [TT] . 

In addition to a-Actinin-4, another known transient 
linker fascin has also been shown to yield network rhe- 
ology consistent with G* oc (iuj) 1 / 2 [18]. In fact, many 
physiological actin cross-linking proteins are dynamic and 
should exhibit a similar cj 1//2 behavior. Indeed, this may 
enable the cell to regulate its response; on timescales 
short compared to the linker unbinding time, the net- 
work is effectively permanently connected — thereby pro- 
viding mechanical resilience — while on longer timescales, 
dynamic linkers allow for complex network flow. This 
ability to flow and remodel is required for many vital cel- 
lular functions, ranging from motility to division. The 
extent to which transient cross-linking affects the me- 



chanical properties of the cell is, however, still largely 
unknown. Interestingly, some rheological measurements 
on living cells have suggested a 1/2 power-law behavior 
on time-scales ranging from several seconds to hours, con- 
sistent with the predictions of our model for a transient 
network [I9j [20]. Further experimental work is needed 
to determine whether this regime is due to the transient 
nature of the cross-links. 
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